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We propose a new technique which enables an event-by-event selection of neutrino-hydrogen inter¬ 
actions in multi-nuclear targets and thereby allows application of hydrogen as targets in experiments 
with neutrino beams without involving cryogenics or high pressure hydrogen gas. This technique 
could significantly improve the reconstruction of the neutrino energy spectra. Since it allows a sep¬ 
aration between hydrogen and the accompanying nuclei, this technique also enables us to measure 
nuclear effects in neutrino interactions directly. 
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The ability to measure the energy spectrum of a neu¬ 
trino [1] beam has many physics implications [2, 3]. The 
accuracy of the measurement depends on the spectral 
shape, the energy reconstruction and the understanding 
of the cross sections of the processes by which the neu¬ 
trino interactions are detected. The conventional mea¬ 
surement of the energy spectra of neutrino beams is via 
charged-current quasi-elastic scattering (CCQE) on nu¬ 
cleons: ly + n ^ l~ + p and D + p ^ + n, where v/v, 

n, p and P stand for neutrino/anti-neutrino, neutron, 
proton and the corresponding charged leptons, respec¬ 
tively. The neutrino energy can be calculated using the 
lepton momentum, assuming a static nucleon in the ini¬ 
tial state [3, 4]. With a nuclear target, the accuracy is 
limited by the binding energy and the Fermi motion (FM) 
of the nucleon, both subject to large fluctuations, and 
by other initial-state uncertainties. When the momen¬ 
tum of the final-state nucleon is measured, the neutrino 
energy can be reconstructed by summing all the final- 
state particle momenta. However, the kinematics of the 
final-state nucleon are altered by final-state interactions 
(FSIs) as the nucleon re-interacts with the cold nuclear 
medium before leaving the target nucleus. FSIs can be 
so strong that the nucleus is excited or even breaks up, 
emitting low momentum particles such as nucleons, pho¬ 
tons and pions, which are stopped near the vertex and 
not detected in tracking detectors, leading to greater bias 
in the reconstructed energy. Therefore a third approach 
(not restricted to CCQF) is to sum the lepton energy and 
the visible energy of the hadronic system [2, 5, 6], which 
is limited by a reduced influence from nucleon initial- 
state uncertainties and by the systematics in measuring 
the energy of neutral particles. Among those initial- 
state uncertainties, multi-nucleon correlations [7-12] are 
under intense study. Such nuclear modifications make 
the calculation of the CCQF cross section difficult. Fur¬ 
thermore, in non-CCQF interactions such as resonance 
production, the final-state pions can re-scatter, exchange 
charge, or be absorbed in the nuclear medium [13]. Such 
background events are often mis-identified as CCQF due 
to the identical final-state particles and therefore intro¬ 


duce an ambiguity in the cross section definition. More 
details about neutrino-nucleus interactions can be found 
in Refs. [14, 15] and references therein. 

Hydrogen is the ideal target for reconstructing the neu¬ 
trino energy because of the absence of these nuclear ef¬ 
fects, however a hydrogen target with high mass is tech¬ 
nically impracticable. In this work, we propose a solution 
which is using a spatial symmetry in the final-state kine¬ 
matics in charged-current (CC) resonance production to 
isolate hydrogen events in targets with a mixture of nu¬ 
clei. This would allow the reconstruction of the beam 
energy spectrum only limited by the knowledge of the 
cross section on hydrogen, which is much better under¬ 
stood than those on nuclei. 

Delta resonances A(1232) can be produced in CC in¬ 
teractions on hydrogen when the neutrino energy is above 
threshold (about 0.34 GeV and 0.49 GeV for Vg and 
respectively). Consider a r'e or interaction on a pro¬ 
ton ly + p l~ + A'''+, where the A++ decays to a 
proton and a positive pion, tt’*'. We define a double- 
transverse axis ztt = Pu X Pi/ \Pi/ x pi\, which is by con¬ 
struction perpendicular to both the neutrino and charged 
lepton momenta, p^, and p). On projecting the proton 
and pion momenta, Pp and p/, onto Ztt, Px-p = Pp ■ ^tt, 
PPT = Ptt'Ztt, one has the double-transverse momentum 
imbalance Sptt = PpT Ptt ^ schematic 

illustration). In the absence of nuclear effects, as is 
expected for a hydrogen target, (jppT is zero, whereas 
i5pTT 7 ^ 0 in the presence of FM and FSI in a nuclear 
target. This is independent of the neutrino energy and 
the resonance kinematics. The ^ppT for a nuclear target 
has the following properties: (1) It is distributed sym¬ 
metrically around zero because the initial proton motion 
and the decay kinematics of the resonance are uncorre¬ 
lated to ipT (except for uncommon cases such as polar¬ 
ized spatially-asymmetric nuclear targets, or if the detec¬ 
tion acceptance varies for different final-state particles). 
(2) Since FM is isotropic, for a given initial proton mo¬ 
mentum Pp (up to about 200 MeV/c for carbon [16], for 
example), (5pTT is broadened from 0 to the same order of 
magnitude as pp. The randomness of FM further smears 
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FIG. 1. Schematic illustration of the double-transverse kine¬ 
matics. The incoming and outgoing particle momenta are 
represented by andp; , Pp andpV, respectively. The double- 
transverse momentum imbalance, 5 ptt is given by -bpxT 
with respect to the axis ztt defined by x p;. 
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FIG. 2. Probability density function (p.d.f.) of the 
double-transverse kinematic imbalance 5 ptt generated by 
NuWro [17] for hydrogen, deuteron, helium, carbon, argon 
and lead targets with neutrino energy 1 GeV. The width of 
the hydrogen distribution is due to the finite bin width. 


out SpTT- (3) The resonance and the decay products ex¬ 
perience FSI. Such modification of the kinematics further 
adds to the broadening of the (ipxT distribution. 

The difference between the (5 ptt shapes for hydrogen 
and nuclear targets is dramatic (see Fig. 2). As can be 
seen in the figure, the shapes of the nuclear distributions 
predicted by NuWro [17] vary only slightly among nuclei 
heavier than deuteron. Therefore for a multi-nuclear tar¬ 
get with hydrogen, assuming perfect detector response, 
one expects a hydrogen signal at (ipxT = 0 on top of a 
symmetric nuclear background that is about 200 MeV/c 
wide. At the reconstruction level, the shape of the hydro¬ 
gen peak, which is still symmetric, is solely determined by 
the detector response. The nuclear background contam¬ 
ination under the hydrogen peak depends on the back¬ 
ground shape and the resonance production cross section 
ratio between the nucleus and hydrogen, which equals 
roughly the atomic number of the nucleus modulo nu¬ 
clear effects. Improving the detector resolution [18] will 
lead to a strong signal enhancement and eventually an 
event-by-event selection of hydrogen interactions. Once 
the hydrogen interactions are selected, the neutrino en¬ 
ergy can be reconstructed by summing the final-state en¬ 
ergy [19]. The energy reconstruction quality is solely de¬ 
termined by detector response and not limited by nuclear 
effects. 

The advantage of JpxT is clear when compared to other 
characteristic variables in the interaction such as the in¬ 
variant mass of A(1232) and the total transverse momen¬ 
tum. The former has an irreducible Breit-Wigner width 
of about 117 MeV/c^ [20] and therefore has no sensitivity 
to reject nuclear background; the latter, which is intrinsi¬ 
cally also zero for a hydrogen target, is asymmetric after 
reconstruction and in general has a long tail resembling 
the Landau distribution [21] due to detector effects. 


For the anti-neutrino interaction on a proton P -|- p ^ 

-I- A°, where A° decays to p -I- tt”, the previous def¬ 
inition and discussions directly apply. This similarity 
enables highly consistent measurements of neutrino and 
anti-neutrino energy spectra. 

Given limited detector resolution, it is important to 
minimize background that has non-zero ^pxT. Because 
FSI can lead to soft nuclear emission, measuring the ver¬ 
tex energy [22, 23] allows tagging and rejecting interac¬ 
tions on other nuclei. And since pion FSI can modify 
the final states, the hydrogen signal purity could be en¬ 
hanced by choosing target materials whose nuclear part 
has a large pion FSI cross section. Non-exclusive back¬ 
ground such as multiple pion production can be rejected 
by vetoing electromagnetic processes and neutral parti¬ 
cles [24]. 

For existing experiments, given a well understood de¬ 
tector response, it may be feasible to perform a combined 
fit to the center region of the JpxT distribution for a 
mixed target, where the hydrogen shape is fixed and the 
background modeling follows the general properties de¬ 
scribed above. If the signal width is at the sub-hundred 
MeV level, the fit may not be sensitive to the detail of the 
complicated nuclear tails. The cross section of the reso¬ 
nance production on hydrogen, which is independent of 
nuclear effects, can be obtained from the signal part. The 
yield ratio between the hydrogen signal and the remain¬ 
ing contribution from the other target nuclei is a pre¬ 
cise measurement of the associated nuclear effects with 
cancellation of detection acceptance and efficiencies for 
both targets. Such measurement of the resonance pro¬ 
duction on hydrogen and nuclei should largely improve 
the understanding of the production mechanism and its 
modification by the nuclear medium. 

More challenging background processes are those with 
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intrinsic zero JpxT- One type of such background is the 
exclusive processes which have identical final states as the 
signal, such as higher mass resonances and non-resonant 
production [25]. Those processes do not affect the neu¬ 
trino energy reconstruction, but make it difficult to define 
the interaction cross section which is needed to deter¬ 
mine the beam energy spectrum. To distinguish among 
underlying processes, the detailed interaction kinematics 
could be used, such as the invariant mass of the hadronic 
system W and the squared four-momentum exchange to 
the proton t. This is feasible thanks to the essential fea¬ 
tures — being exclusive and nuclear effect-independent 
— that make the neutrino energy reconstruction precise. 
Because the cross sections have different dependence on 
W and t, efficient separation should be possible. Since 
the relevant final state kinematics do not depend on the 
identity of the intermediate state, alternatively one could 
extend the definition of the production channel to include 
all contributions that have exclusive p7r+ final states and 
calculate the corresponding cross section [26, 27]. An¬ 
other type of background is the contamination 

in the (z/^) CC interaction. Such “wrong sign” back¬ 
ground occurs more often in an anti-neutrino beam cre¬ 
ated in a proton accelerator. Because of the similar par¬ 
ticle identification (PID) signals for the muon and pion, 
the fi~ +TT^ +i> final states from the background might 
be misidentified in the event selection as 7r“ -|-/r+-l-p from 
P^, and vice versa. For a large enough detection volume, 
the properties of stopped pions can be used to enhance 
the pion identification. In addition to measuring trajec¬ 
tory “kinks” [28] and the Michel electrons [29], negative 
pions that are stopped and absorbed by nuclei give rise to 
soft nuclear emission which can be measured by calorime¬ 
try [30]. The lepton, 7r+ and tt~ signatures are different 
and serve as an important tool not only for 1/tt separa¬ 
tion, but also for rejecting “wrong sign” contamination. 
In addition, variables like W, which is calculated from 
the true p/i system due to wrong PID, can provide back¬ 
ground rejection power due to the unphysical kinematic 
combination [31]. 

In this work, MC simulation of the T2K ND280 detec¬ 
tor [32] is used to demonstrate the measurement of CC 
resonance production in neutrino-hydrogen interactions, 
zz^ -I- H —+ A++. Neutrino interactions in ND280 
are simulated using the Neut event generator [33]. The 
event reconstruction uses the first fine-grained detector 
(FGD [34]) as a plastic scintillator (polystyrene) inter¬ 
action target and the neighboring gaseous time projec¬ 
tion chamber (TPC [35]) downstream to measure the mo¬ 
menta and specific energy loss {dE /dx) of the final-state 
particles. There are about 300 true A++ events on hy¬ 
drogen in the neutrino energy range 0.5 < < 5 GeV 

in the reconstructed sample. 

The distribution of the reconstructed double- 
transverse kinematic imbalance 5ptt in the ND280 
acceptance (detection efficiencies apply) is shown in 
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FIG. 3. Distributions of the reconstructed (JpxT in the ND280 
acceptance (detection efficiencies apply) simulated for hydro¬ 
gen and carbon target nuclei in the neutrino energy range 
0.79 - 1.26 GeV. Vertical bars are MC statistical errors. The 
hydrogen signal is fit to a Cauchy distribution. The carbon 
distribution is area-normalized to the hydrogen. 


Fig. 3. Due to the absence of nuclear effects, for hydrogen 
the resolution is determined by the detector response, 
which can be described in the simulation by a Cauchy 
function 1/N ■ dN/dSpTT = I/tt • + (<5ptt — w)^] 

with the event count A, mean m and width tr. The 
resolution is shown to be significantly smaller than the 
nuclear broadening in the carbon target nuclei. The 
nuclear rejection factor defined as the efficiency ratio 
between hydrogen and carbon is about 3.7 (2.3) in the 1 
(3) a interval. The reconstruction performance of ^pxT 
as a function of the neutrino energy is shown in Fig. 4. 
It indicates a better hydrogen selection at lower energy. 

From the MC true hydrogen-resonance events in the 
reconstructed sample, the neutrino energy E^^ is calcu¬ 
lated directly using the kinematics of all final-state parti¬ 
cles. The uncertainty in an example Ey^ bin is shown in 
Fig. 5. Like (5pxT, it is determined by pure detector res¬ 
olution due to the free and static proton target and can 
be described by a Cauchy function. As a comparison, 
the neutrino energy in the same range is reconstructed 
with all final-state kinematics from CCQE interactions 
on carbon assuming static neutron targets. This simple 
reconstruction suffers from nuclear effects in addition to 
the detector resolution and is biased by about —2% with 
a large spread. In addition, the neutrino energy distribu¬ 
tion reconstructed only with the muon kinematics [4, 36] 
shows worse performance. The comparison between both 
CCQE methods indicates that, for the relatively weak 
FSI predicted in the current model, an independent mea¬ 
surement of the proton kinematics help improve the Ey 
resolution even in the presence of FSI. The simulated Ey^ 
detector response (the Cauchy mean and width) is fur¬ 
ther shown as a function of the true energy in Fig. 6. The 
energy scale is seen to be constant with a bias of about - 
























4 



FIG. 4. The Cauchy mean (upper) and width (lower) of the 
reconstructed 5ptt as a function of the neutrino energy. Ver¬ 
tical bars are MC statistical errors, while horizontal error bars 
stand for the bin span of the neutrino energy. A fit of a con¬ 
stant is applied to the mean values. 


1%, and the resolution slowly increases from about 3% to 
10% in the 0.5 - 5 GeV region. The worsening of resolu¬ 
tion at high energy is general for a tracker measurement, 
in contrast to calorimetry [2, 5, 6]. 

Given the fact that ND280 was optimized to measure 
GGQE interactions and designed to achieve a perfor¬ 
mance that was to be limited by the nucleon initial-state 
uncertainties in the nuclear target [34, 35], alternative 
optimization and state-of-the-art technology may allow 
better performance. Since the GG resonance production 
cross section rises to a maximum at a neutrino energy 
greater than about 3 GeV [20], the high energy neutri¬ 
nos produced in the NuMI [37] and LBNF [2] beam lines 
are optimal to realize the proposed method. It would 
be challenging, yet very attractive, to combine the pro¬ 
posed use of hydrogen targets and the liquid argon TPC 
projects [2, 38-41] which are to realize superb tracking 
and calorimetry performance on a massive scale. Fi¬ 
nally, it would be interesting to demonstrate the method 
with anti-neutrino beams as well as for Vslve interac¬ 
tions in oscillation appearance and investigate the impact 
on the physics programs in future experiments, consid¬ 
ering also the advantage of the identical target, interac¬ 
tion kinematics and phase space in the proposed channels 
vjv -b p —4- A -T -I- p -I- TT^, in comparison to the 
conventional GGQE interactions. 

The authors would like to express their gratitude to the 
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FIG. 5. Deviation of the reconstructed neutrino energy from 
the true value. Reconstruction with the CG resonance pro¬ 
duction on hydrogen targets is compared to the ones using 
GGQE on carbon: one summing both muon and proton re¬ 
constructed kinematics (“p -bp”), the other using the muon 
to calculate the proton momentum assuming a static initial 
neutron (“p + ”) [4] . The carbon distributions are area- 

normalized to the hydrogen. 



FIG. 6. Neutrino energy scale (upper) and reconstruction 
resolution (lower) via GC resonance production on hydrogen. 
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